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Abstract 

Q^l The Gell-Mann and Low switching allows to transform eigenstates of an unperturbed 

Hamiltonian Hq into eigenstates of the modified Hamiltonian Hq + V. This switching can 
be performed when the initial eigenstate is not degenerate, under some gap conditions 
■ with the remainder of the spectrum. We show here how to extend this approach to the 

case when the ground state of the unperturbed Hamiltonian is degenerate. More precisely, 
we prove that the switching procedure can still be performed when the initial states are 
eigenstates of the finite rank self-adjoint operator VqVVo, where Po is the projection onto 
1^ . a degenerate eigenspace of Hq . 

m 

oo ■ 1 Introduction 

^ ■ Adiabatic switching is a crucial ingredient of many-body theory. It provides a way to express 

^ . the eigenstates of a Hamiltonian Hq + V in terms of the eigenstates of Hq. Its basic idea 

Q I is to switch very slowly the interaction V, i. e. to transform Hq + V into a time-dependent 

Hamiltonian of the typical form Hq + e"^'*'!/, where the small parameter e > eventually 
vanishes. It may be expected that an eigenstate of Hq + V is obtained by taking the limit 



! of an eigenstate of Hq, evolved according to the time-dependent Hamiltonian Hq + e~^^^W 

I when e tends to zero. It turns out that this naive expectation is not justified since the time- 

dependent eigenstate has no limit when e — > because of some non-convergent phase factor. 
When the initial state belongs to a non degenerate eigenspace, Gell-Mann and Low solved the 
problem by dividing out the oscillations by a suitable factor [7]. The ratio becomes, in the 
limit e ^ 0, the Gell-Mann and Low wavefunction. Mathematically, the convergence of this 
procedure has been proved in 1989 by Nenciu and Rasche [16], elaborating on the adiabatic 
theorem [3l[T2l[6]. 

On the other side, the physics community realized about fifty years ago [2] that a general- 
ization of the Gell-Mann and Low formula is needed in the case of a degenerate eigenvalue of 
Hq. This happens in many practical situations, for instance when the system contains unfilled 
shells. This problem has been discussed in several fields, including nuclear physics, solid state 
physics, quantum chemistry and atomic physics, see the references in [11[5]. In most cases, it 
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is assumed that there is some eigenstate in the degenerate eigenspace Sq of Hq for which the 
Gell-Mann and Low formula holds. In general however, the Gell-Mann and Low formula is 
not applicable when this state is chosen at random in the degenerate subspace, as illustrated 
in the simple model analytically studied in [4]. 

We show in this paper that the switching can be performed provided the initial eigen- 
states are also eigenstates of VoWo\g^^, the perturbation restricted to act on the degenerate 
eigenspace. If the latter operator has itself degenerate eigenvalues, a further analysis is re- 
quired, as discussed in Section 13.41 The result is based on the recent progress in the math- 
ematical analysis of adiabatic problems, see [151 E E Ell El [HI [13 [22l [10] and references 
therein. 

The physical consequences of our result are discussed in the companion physics paper [5], 
where we also comment on the formal relation with different types of Green functions. 

Acknowledgements: G.P. is grateful to S. Teufel and J. Wachsmuth for a useful discussion 
in a preliminary stage of this work. We gratefully thank an anonymous referee for useful 
comments and remarks, which encouraged us to generalize the result appearing in the first 
version of the paper. 



2 Statement of the results 

2.1 Spectral structure of the problem 

Consider a Hilbert space 7i, a self-adjoint operator Hq, with dense domain D{Hq) C TL, and 
a symmetric perturbation V , ffo-bounded with relative bound a < 1. Then, according to 
the Kato-Rellich theorem (Theorem X.12 in [E]), Hq + W is self-adjoint on D{Hq) for any 
< A < 1. We denotes 

H{X) = Ho + XV, 

with A € [0, 1]. In all this study, we assume that the spectrum has the following structure. 

Assumption 1 (Structure of the spectrum). The spectrum of H{X) = Hq + XV, X € [0,1], 
consists of two disconnected pieces 



a{H{X)) = aN{X) U (a(if(A))VAr(A) ^ 
where cr7v(A) is a finite subset of the discrete spectrum: 

aNiX) = {Ej{X), i = 1, . . . , at} c adisc (^'(A)) , 

and the initial state is degenerate: Ej{0) = Ek{0) for all I < j, k < N . 

In order to apply results and techniques from adiabatic theory [31 [121 [111 [I] , we make the 
following standard assumption on the existence of a gap in the spectrum. 

Assumption 2 (Gap condition). There is a gap between the two parts of the spectrum, in the 
sense that: 

A(A) = _inin^(^min| Ej{X) - E , E e a{H{X))\{Ei{X), . . . ,En{X)]^ 



^ For reasons that will become clear once a time variable is introduced, we will always denote with a 
functions of the variable A e [0, 1]. Untilded functions will have time as an argument. 
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is bounded from below by a positive constant: 

inf A(A) = A* > 0. 

AG[0,1] 

The projectors associated with the N eigenvalues Ej{\) (counted with their multiphcities) 
are denoted by Pj{X), for 1 < j < M with M < N. The projector onto the subspace 
orthogonal to the eigenspace spanned by the N eigenvectors is Pn+i{X) = I — X^jli-Pj(A). 
We denote in the sequel 

M 

the projector onto the eigenspace £o = Ran(7'o) spanned by the N degenerate eigenstates of 
Hq. For simplicity, we assume that the perturbation V is sufficient to split the degeneracy (so 
that M = A^), in the sense that the following assumption holds true. 

Assumption 3 (Degeneracy splitting). The finite rank self-adjoint operator VqVVq : Sq ^ Sq 
has non-degenerate eigenvalues, and there is a gap between the N first levels in the interval 
(0, 1]; for any X* > 0, there exists a (depending on X*) such that 



inf min 

A*<A<1 k^l 



Ek{X) - EiiX) 



>a>0. (2.1) 



This implies that the projectors -Pj (A) are rank-1 projectors for any A > (since it can be 
proved that the perturbation V is enough to split the eigensubspaces, and the gap condition 
on (0, 1] ensures that no crossing can happen; see Section [STTl for more details). 

Remark 4. Assumption\^may be relaxed in several ways. First, the operator VqVVq can have 
degenerate eigenvalues, but then higher order terms should be considered in the perturbative 
expansion of the eigenvalues. The gap assumption can be relaxed as well, and some crossings 
could be allowed. Besides, the general case of M < N projectors of ranks greater or equal 
to 1 can be treated similarly upon modifying the condition PjiX) ~ Pj{^) < 1 required in 
Theorem\7\ below. All these extensions are discussed in Section [X^ 

2.2 Switching procedure 

Consider, for r G ( — oo,0], 

H{T)=H{f{T))=Ho + f{T)V, 

where the switching function / has values in [0, 1] (in order for the operator H{t) to be 
well-defined as a self-adjoint operator on D{Ho)). We denote by Pjir) the eigenprojectors 
and eigenvalues corresponding to the first eigenvalues Ej{T) of H{t); also, i-*Ar+i(r) = 
^-ELiPkir). Of course, 

Pj{t) = P,{f{T)), Ej{t) = E,{f{T)). 
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For the subsequent analysis, we assume that 



Assumption 5. The switching function f : (— oo,0] —^ [0, 1] is a function such that 
(i) f is non- decreasing; 
(ii) /(O) = 1 and hm /(r) = 0; 

r— oo 

(ill) /,/"€LH(-oo,0]). 

The most common choice in practice is /(r) = e"^. Notice however that any non- 
decreasing compactly supported function with /(O) = 1 satisfies the above assumptions. In 
the latter case, the monotonicity of / implies that the support of / is a compact interval 
[i?j-,0], and f{t) > for t € (i?j,0]. The assumption / G ensures that the adiabatic 
evolution (see (|3.13|) below) is well-defined. 

As a consequence of these assumptions, /' > and /' € L^((— oo, 0]) nL°°((— oo, 0]), hence 
/' G L^((— 00,0]). Indeed, the boundedness of /' is a consequence of the fundamental theorem 
of calculus and the fact that /" G l^{{-oo, 0]). Besides, /° /' = /(O) - f{t) < 1, and /' > 0, 
hence /' G L^((-oo,0]). 

Remark 6. It can be shown that eigenprojectors and eigenvectors are analytic with respect to 
A = /(r) (see Section \3. When the switching function f is analytic, the eigenvalues Ej{T) 
(and the associated eigenvectors and eigenprojectors) are also analytic with respect to t. 

We denote by Us{s, sq) the unitary evolution generated by H{es), i. e. the unique solution 
(which is well-defined by Theorem X.70 in |18j l of the problem: 

i^^^^^j^ = H{es)U,{s,so), Us{so,so)=l. 

In order to remove divergent phase factors (see the proof in Section [3.3. it is convenient to 
consider evolution operators in the interaction picture: 

f/e,int(s,5o) = e'^^°C/,(s,so)e-^^°^°. 

It is actually more convenient to rescale the time and to consider a macroscopic time t = es. 
The unitary evolution U^{t,tQ) in terms of the macroscopic time is the solution of 

ie^^^^^^ = H{t)U%t,to), U%to,to)=I, 

and, in the interaction picture, 

U[^,{t, to) = e'*^°/^ U%t, to) e-'*"^o/^ 

Standard results show that U-^^^{t, —oo)^p = lim(„_+_oo f^hit(^' *o)^ exists for tp G D[Ho) (for 
instance, by using a Cook's type argument and rewriting this operator as the integral of its 
derivative with respect to to)- 
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2.3 Main results 



We are now in position to state our main results. 

Theorem 7. Suppose that the gap conditions on Hq and V (Assumptions [I\ andl^ are sat- 
isfied, and that the perturbation term V lifts the degeneracy (Assumption Consider a 
switching function verifying Assumptions^ Let (■01, • • • , V'iv) an orthonormal basis of Sq 
which diagonalizes the bounded operator VqV'Pq\^^. Then, if 

||P^.(_oo)-P,(0)||<l, (2.2) 

the limit 

^, = lim , ^^t(O, -oo)0, 
' .^0(0,1 C/^^(O,-oo)0,) 

exists and is an eigenstate of Hq + V . 

Notice that, for a generic state € RanPo which is not an eigenvector of 'Po^'Pol^-^ the 
above hmit generically does not exist, as showed in by using a simple toy model. It is 
therefore crucial to select the appropriate initial states, so that the Gell-Mann & Low limit 
(12.31) does exist. 



As an intermediate step, the eigenprojector Pj(0) and a corresponding eigenfunction 
can be recovered by Kato's geometric evolution [T2] . 

Definition 8. The Kato evolution operator A[s,sq), for s, sq £ K is the unique solution of 
the problem 

dA{s,SQ) 
ds 

with 



K{s)A{s,sq), A(so,so) = I, (2.4) 



K{s) = -Y^P,{sf-^{s). 

By our assumptions, the operator K[s) is uniformly bounded (see the comment after Defini- 
tion [TT]). The Kato evolution operator is a unitary operator which intertwines the spectral 
subspaces of H[s) and H{so), in the sense that 

A{s,so)Pj{so) = Pjis)A{s,so). 

Equipped with this notation, we have the following result, where no condition analogous 
to (12.21) is assumed. 



Proposition 9. Let Assumptions{^ and\^ be satisfied. Let {ipi, . . . , tpj\f) be an orthonor- 
mal basis of So which diagonalizes the operator 'PqV'Pq\^^. Then 

:= A{0, -oo)V'i 

is an eigenvector of Hq + V. 
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It is actually much simpler to consider the geometric evolution operator A rather than the 
evolution operator C/^^ since less conditions on Hq and V are required. Indeed, there is no 
denominator which needs to be considered in order to remove a divergent phase. However, 
the many-body theory used in physics is defined in terms of U[^^^ and not in terms of A. 

We sketch shortly the structure of the proof, which is done in four steps: 

(i) first, we use the Kato geometric evolution backward in time, in order to identify, though 
in a non explicit manner, the initial subspaces of Vq whose vectors can be considered as 
convenient initial states; 

(ii) in a second step (Section l3.2p . we give an explicit description of these initial subspaces, 
in terms of the eigenvectors of ^o^^ol^;^- At this stage, we are already in position to 
prove Proposition [9l 

(iii) then, we show how the limit of the full evolution ^7^^. can be related to the geometric 
evolution as e — > (Section [33]) • A first step is to introduce an intermediate concept, the 
adiabatic evolution, which takes some dynamics into account (arising from the Hamilto- 
nian operator). The adiabatic evolution is also an intertwiner. Since intertwiners differ 
only by a phase (in a sense to be made precise), and, provided this phase can be removed, 
the adiabatic evolution can be reduced to the geometric one (see Section [3.3.1|) : 

(iv) the last point is to show that the limit as e — > of the full evolution is the adiabatic 
evolution (see Section [3.3.2p . 

Steps (iii) and (iv) are straightforward extensions of previous results in adiabatic theory, and 
we heavily relied on the paper by Nenciu and Rasche [16] for Section 13.3.11 and the book by 
Teufel [22] for Section [3X2l 

3 Proof of tlie results 

3.1 Geometric evolution and definition of the initial states 

In view of the local gap assumption, the projectors and eigenvalues of H{X) are real analytic 
functions of A G (0, 1]. Besides, Theorem II. 6.1 in [13] shows that the eigenvalues Ej and pro- 
jectors Pj can be analytically continued in the limit A ^ 0. The Kato construction of unitary 
operators A intertwining projectors can then be performed, see for instance Theorem XII. 12 
in [19] or Sections II. 4 and II. 6. 2 in [13]. Consider the operator 



Since K{X) is uniformly bounded, thejaperator A{X, Aq) is well-defined and strongly continuous 
(see Theorem X.69 in [IB]). Besides, ^(A, Aq) is unitary, and intertwines the spectral subspaces: 



K{X) 




first proposed in [12], and the unique solution of 



dA{X,Xo) 
dX 



K{X)A{X,Xo), A{Xo,Xo)=l. 



(3.1) 



P,(A)=^(A,Ao)P,(Ao)^(A,Ao)*. 
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It is also easily shown that A{X2, Xi)A{Xi, Xq) = A{X2,Xo), for instance by computing the 
derivative of both expressions with respect to Ai and using the uniqueness of the solution 
of (IXTD . 

We define the initial subspaces by evolving backwards eigenstates of the Hamiltonian H{X) 
for which the perturbation has split the degeneracy: the corresponding eigenprojector is de- 
fined as 

P;-'':=AiO,X)P,iX)A{X,0), (3.2) 

the definition being independent of A > 0. 

Eigenstates of H{1) = Hq + V are then obtained by evolving initial states belonging 
to the range of Pj"^'* according to the Kato evolution operator. Indeed, j4(l,0)Pj"* = 

A{l,0)A{0,X)Pj{X) = A{l,X)Pj{X). Thanks to the intertwining property of A, it holds 

Pj{l) = I(l,0)Pf '*I(0, 1). (3.3) 



3.2 Characterization of the initial states 

The above paragraph shows that it is crucial to identify Ran(Pj°'*). We now characterize 
these spaces by an explicit condition. 



General expressions of the eigenvalues and eigenvectors. Since the eigenvalues and 
eigenprojectors of H{X) are analytic in A € [0,1], the following expansions are valid for 1 < 
j<N: 



+ 00 



E,(A) = 5]A"i?,,. 



(3.4) 



and 



n=0 



+ 00 



n=Q 



Of course, Ej^ = Eq = Ej{0), the common value of the energy in the degenerate ground-state. 
Notice also that the operators Pj^n are not necessarily orthogonal projectors. 

To define Pj{X), it is more convenient to consider an eigenvector 4>j{X) associated with 
Ej{X), i. e. a non-zero element of Ti. satisfying 



H{X)cl,,{X) = E,{X)<P,{X). 



(3.5) 



Such an eigenvector can be chosen to be analytic, by the same results which allow to conclude 
to the analyticity of the eigenprojectors. We therefore write 



(3.6) 



n=0 



Once such an eigenvector is known, the analytic eigenprojector can be constructed as 



wMm I \ um\ 
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The aim of this section is to provide an exphcit expression of the leading terms of the 
above expansions, in order to have a more exphcit definition of To this end, we first 

construct a basis of Sq, which will turn out to be particularly useful to characterize the terms 
in the expansions (|3.4p and (13. 6p . 

Diagonalization of VqVVq. Since VqWo and Vq commute, it is possible to construct an 
orthonormal basis {ipifi, ■ ■ ■ , ^n,o) of Sq such that 

VoWo ^j,o = ctjipjfi (3.7) 

for some real numbers ai, . . . , a^r, and 

Vj ^ k, { ifkfi, VoVVo ifjfi ) = 0. (3.8) 

Expressions for the terms in the expansions (|3.4|) - (|3.6|) at order 1. We identify the 
terms associated with the same powers of A in (13. 5j) . An additional normalization condition 
should be added in order to uniquely define the solution, so we impose 

VAg[0,1], {ip,fl,c^j{\)) = l, (3.9) 

as is done in [20]. As will be seen below, this condition is simpler to work with than the 
standard condition ||(/)j(A)|| = 1. The identification of the terms in (|3.5p gives, for 1 < j < N , 
the following hierarchy of equations: 



(^0 


- Eo)ipjfl 


= 0, 




{Ho 


- Eo)(pj,i 


= {Ej,i 




{Ho 


- Eo)(pj,2 


= {Ej,i 


- V)(fj^i + Ej^2'Pj,o 



and, for n > 2, 

n-l 

{Ho - Eo)ipj^n+l = {Ej^l - V)ipj^n + ^ Ej^n+l-m^j,m- (3.10) 

m=0 

The equation on the term of order zero does not give any information on the choice of the 
initial states ipj,o- This information can be obtained from the first order condition: 

{Ho - Eo)ipj,i = (^,-,1 - V)ip,fl. (3.11) 

A necessary condition for this equation to have a solution is that the right-hand side belongs 
to £^ (since the left-hand side does): 

VI <j,k<N, { ^k,o, (%i - V)ipj,o ) = 0. (3.12) 

This requires 
and 

VA:/i, {ipk,o,Vipjfl) = 0. 
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Therefore, the conditions (|3.12j) for k ^ j cannot be fulfilled for a general basis. A necessary 
condition is that the basis {'^k,o}k=i,...,N of fo diagonalizes VqVVq. Besides, the first-order 
term in the energy shifts are exactly the eigenvalues of VqVVo- This condition determines 
uniquely the basis when VqWo has non-degenerate eigenvalues. If this is not the case, infor- 
mation about the higher order equations in the hierarchy is needed (see Section [37 



Remark 10. Assuming that the bands do not recross after the initial splitting, and if the 
degenerate state is the ground state of Hq, then the ground state of Hq + V is obtained by 
following the eigenstate associated with the lowest Ej^i. 

Once the initial basis and the first energy shifts have been defined, the first order term 
in the variation of the eigenstates can be obtained from (|3.11|) as the sum of the reduced 
resolvent applied to the right hand side, and some solution of the homogeneous equation 
{Ho - Eo)i' = 0: 

N 

'^i.i = ^],k'Pk,o + (Ho - Eq)-^ {Ej^i - V)ipj,o 
k=l ^0 



where 



Ro = {Ho - Eor^ ^ , = (I - Vo) {Ho - Eo)'^ (I - Vo) 



is a bounded operator from £^ to n D{Ho), and cj^ = in view of the normalization 
condition (|3.9|) . The coefficients c\ j (for k ^ j) are undetermined at this stage. They have to 
be chosen so that the right hand side of the next equation in the hierarchy is in £f^. 

Conclusion: characterization of the initial subspaces. The above computations show 
that Pj{\) = Pjfl + 0(A), with Pj^ = \ipjfi){ipjfi\. Besides, P(0,A) - I|| = 0(A) in view of 
the differential equation (|3.ip satisfied by A. The initial subspace (13. 2p is therefore 



A{0, \)P,{\) = hm ^(0, A) [P,- + 0(A)] = P,- q. 

A — ^0 



Proof of Proposition [9l Let ip & So an eigenvector of VoVVo- Then, there exists 
j e{l,...,N} such that ip G Ran(P,-o) = Ran(P^'^i*). Using ([33]), it follows 



^(0, -oo)^j = A{1, 0)'ipj e Ran (^Pj(l)j , 

which proves the claim. 

3.3 Adiabatic evolution and limit of the full evolution 

Definition 11. The adiabatic evolution operator Ua{s,so) is defined for (s,so) € as the 
unique solution of the problem 

.dUA^s^ = Ha{s)Ua{s, so), Ua{so, so) = I, (3.13) 
as 
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where the adiabatic Hamiltonian is 

HA{s) = H{s) + iK{s), 

with 

Notice that K{s) = f'{s) K{f{s)) so that 

\\K{s)\\<Cf'{s) (3.14) 

for some constant C > 0. Therefore, K{s) is uniformly bounded since /' is bounded by our 
assumptions on the switching function. 

Compared to the geometric evolution (13. Ih . a Hamiltonian term has been added, which 
is at the origin of some dynamical phase factor in the dynamics. The adiabatic dynamics is 
well-defined in view of the assumptions made on Hq, V and / (see Theorem X.70 in [IHj). A 
simple computation shows that it intertwines the spectral subspaces: 

Pj{s) = UAis,so)Pj{so)UAis,soy. 

Switching to the interaction picture, we define 

f/A,int(s, so) = e'^^° UAis, so) e-^^«^o. 

The factor e is introduced by slowing down the switching as 

.dUe,Ais,So) 



ds 



HA{es)Ue,A{s,so), Ue,A{so,so) =1, (3.15) 



and the corresponding operator in the interaction picture is e'*^°C/e^A('S, sq) e '^o^o. It is 
convenient to rewrite the evolution (|3.15l) in the rescaled time variable t = es: 

ie^^^^^ = Hi{t)Ui{t,to), Ui{to,to)=I, (3.16) 

with Hj^{t) = H{t) + ieK{t). The associated operator in the interaction picture is 

UUntit,to) = ei*^«/^t/i(t,to)e-"°^°/^ 
Theorem [7] is then a consequence of the following results. 

Lemma 12. Let ipj G pj^it ^dg^ji^ed by (|3.2I) ). Then, under the assumptions of Theorem\% 
the vector 

f^A,int(0'-Oo)^i C/A,int(0,-Oo)^j 



{^3 I f^A,int(0' (V'j I f^A,int(0, -Oo)i^j) 

is an eigenstate of Ho- 

Lemma 13. Let ipj S P™*. Then, under the assumptions of Theorem\7[ 



(3.17) 



C/^,(0,-oo)^,- Ul.^,{0,-oo)i^j 



lim "intv"^ -^yj ^'"'"^ '—^ = 0. 

V ( V-, I ^fnt(0, -OO)^',) ii^j I UlUO, -00)i^j) 
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3.3.1 Proof of Lemma [l2] 



We show first in tliis section that ^pj can be transformed into an eigenstate of H(0) = H{1) 
using the adiabatic evolution defined from (I3.13j) , and then the equahty of the ratios (|3.17l) . 
The proof presented here reproduces the argument of Nenciu and Rasche [16], which was 
given in the case = 1 with our notation, but can be appHed mutatis mutandis to the case 
considered here. We however present the proof for completeness. 

Evolution in the case e = 1. Since both Ua and A are intertwiners, they differ only by a 
phase which commutes with the spectral projectors. Indeed, define 

So) = ^(s, so)*Ua{s, So), 
so that Ua{s, So) = A{s, sq) <1>(s, sq). Then, 

[$(s,so),P,(so)] =0, 

as can be seen using the intertwining properties: 

[$(s, so), P;(so)] = A{s, so)*Ua{s, so)Pj{so) - Pj{so)A{s, so)*Ua{s, sq) 

= A{s,soTPj{s)Ua{s,so)-A{s,so)*Pj{s)Ua{s,so) = 0. 



The time-evolution of the phase matrix can be simplified due to this commutation property. 
First, 

d$(s,so) 



ds 

since K{s)* = —K{s). Besides, 



-iA{s,soTH{s)UA{s,so), 



/N+l \ /N+l \ N+l 

Hs,so) = Pk{so) $(s,so) ^-^(^0) = E ^k{s,so), 

\k=l J \k=l J k=l 

where <I>fc(s,so) = Pkiso)^is, so)Pk{so). The time evolution of the projected phase-matrix is 
a scalar phase since 

— $fc(s,so) = -iEk{s)^kis,so), 

hence 

$(s,so)Pj(so) = exp ^-i^ Ej{r)dr^ Pj{so). 

The geometric evolution and the adiabatic evolution are therefore related through some global 
dynamical phase: 

[/a(s,so)P,(so) = ^(s,so)^>(s,so)Pj(so) = exp ^-i^ Ej{r)dr^ A{s, so)Pj{so). 



To describe the asymptotic evolution, we follow closely the approach of [16]. In order for 
Ua{s, SQ)Pj{so) to be defined in the limit sq — oo, it is important to work in the interaction 
picture. Then, 

UA,int{s,so)Pj{-oc) = e'^^M(s,so)$(s,so)e-'^o^op,(-oo) 

= e-''°^'>e''^°A{s,so)e-''^°e''^°^{s,so)Pj{-oo). 
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Using 

^{s,so)Pj{so) = Pj{so)'^{s,so)Pj{so) = exp ^-i^ Ej{r)dr^ Pjiso), 

it holds 

= exp (^-i^' E,{r) dr - \soE^ e'^^'^Pjiso) + e-'^°^Oe^'^^°$(s, so)(Pj(-oo) - Pj{so)) 

= exp (^-i^' E,{r) dr - iso^o) [e*^^"Pi(-oo) + e'^^°(Pj(so) - ^,(-00))] 

+e-''°^'^e''''°^s,so){Pji-oc) - P,{so)) 
= exp(^-i^ Ej{r) - Eodr^ P,{- 00) + W{s, so) {Pj (so) - Pj{- 00)), 

where \\W\\ < 2. Since A ^ Ej{X) is on the compact interval [0, 1], there exists a constant 
C > such that 



|£;,-(r) 



Ej{f{r))-E,{0) <Cf{r). 



Since / € L^((— 00, 0]), this shows that the function r i-^- Ej{r) — Eq is integrable on (—00, 0]. 
Besides, P{so) — > 00) when sq — > —00. The limit sq — > —00 of ?7A,int(S) •so)-Pj(— cjo) is 
therefore well-defined: 

UA;mt{s, -oo)Pj(-oo) = exp ^-i j Ej{r) - E^dr^ e''-^M(s, -oo)e-^'^«Pj(-oo). (3.18) 
The above equality reads, for s = 0, 

t/A,int(0,-cx))Pj(-oo) =exp(^-i j Ej{r) - Eodr^ A{0,- 00) Pj{- 00). 

Since P,(0)^(0, -00) = A{0, -oo)P,(-oo), it holds, for tpj £ Pj"* = Pj(-oo) = R&n{ipj^o), 

Pj{0)i'j = ^(0, -oo)P,- (-00)^(0, -oo)>j = A{0, - 00)* il^j) ^(0, -oo)^j. (3.19) 
Finally, 

^,■(0)^^^ ^ Pj(0)^i ^ ^(0, -oo)iPj ^ UA,int{0,-oo)i'j 

\\Pj{0)lPj ||2 ( I P,(0)^j ) ( I A(0, -00)^j ) ( V', I ?7A,int(0, -oo)Vj ) ' 

which shows that the adiabatic evolution transforms the initial eigenstate into an eigenstate 
of H{1) provided ||Pj(0)V'j|| 7^ 0, which is the case when \\Pj{0) - Pj(-oo)|| < 1. 

Evolution in the case e > 0. Let us conclude this section by proving the equality (|3.17l) . 
Computations similar to the ones performed in the case e = 1 lead to 

C/i,i,t(0, -oo)Pj(-oo) = exp (^-'-j Ej{T) - Eq dr^ A{0, -oo)Pj{-oo). 

This can be seen for instance by noticing that (|3.16|) can be rewritten in the form (|3.13p . upon 
considering the Hamiltonian H/e in (|3.13j) . Therefore, U]^^^^{0,— 00) Pj{— 00) is equal, up to 
the e-dependence in the phase factor, to f/A,mt(0, — oo)Pj(— 00). The non convergent phase 
factor can be eliminated precisely by considering the Gell-Mann and Low ratio (|3.17|) . 
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3.3.2 Proof of Lemma [l3] 



It is sufficient to prove that 

lim ||;7^(0, -oo) - ;7i(0, -oo)\\ = 0, 

which indeed gives the result since 

\\U[^,{t,to)-Ul,^,{t,to)\\ = \\U'{t,to)-Ui{t,to)\\. 

Notice that, although none of the operators C/^(0, — oo), Uj^{0, —oo) has a limit when e — > 0, 
the difference goes to in this limit. 

The proof is based on the proofs of Theorem 2.2 and Corollary 2.5 in the book by 
Teufel [22], which are extended to the case of non-compactly supported switching functions 
and > 1 with our notation. In this section, C and C denote constants, which may change 
from line to line, but are always independent of t, e, etc, and depends only on the relative 
-ffo-bound of V, on A^, on A* and on bounds on the functions Pj and their derivatives on 
[0,1]. 

We denote by Sj{t) > the local gap around Ej{t): 

5j{t) = mm i^\Ej{t) - E\, E £ a{H{t))\{Ej{t)}y 

Notice that Sj{t) > when f{t) > 0, but 6j{t) when f{t) — > since the initial eigenvalue 
is A^-fold degenerate (see Assumption [3]) . In fact, the analysis of Section [3?2] shows that there 
exist ai , 02 > such that 

6,{t) 



Oil < 
when f{t) > 0. 



fit) 



< aa (3.20) 



Rewriting the difference as an integral. The difference between the two unitary evolu- 
tions is rewritten as the integral of the derivative, as: 

U%t,to)-Uiit,to) = -U%t,to) f ^{U%to,t')Ul{t',to)) dt' 

J to 



--U'{t,to) fu'{to,t') [H{t')-HA{t')] Ui{t',to)dt' 

£ Jto 

-U%t,to) [ U'{to,t')K{t')Ui{t',to)dt'. 



Jto 

The idea is to rewrite K{t) as a commutator, so that t i— U^(to,t)K(t)U]^{t, to) is the derivative 
of a function (up to negligible terms), and an integration by parts gives the required estimates. 
The proof proposed here is an extension of the proof presented in [22], Chapter 2] in the case 
when several pieces of the discrete spectrum are considered independently. It would also have 
been possible to use the twiddle operation introduced in [T], which is, in some sense, the 
inverse operation of the commutator with the Hamiltonian. 
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Construction of the function used in the commutator. Consider — oo < t < such 
that f{t) > (for compactly supported switching functions, this means that t is in the interior 
of the support). Define 

1 \ 



with, for 1 < j < iV 

2i7r 
1 
2br 



Fjit) = ^ (f) PHt)R{z,t)R{z,t)dz, (3.21) 



GAt) = — d) R{z,t)R{z,t)P^{t)dz, (3.22) 



where 



R{z,t) = {H{t)-z)-\ Riz,t) = I [iH{t)-z)-'] = -R{z,t)^^Riz,t), 

and Tj{t) is a contour enclosing Ej{t) and no other element of the spectrum (such a contour 
exists in view of Assumption [3|) . For j = + 1, we denote by r7v+i(t) a contour enclosing all 
the first N eigenvalues Ek{t), k = 1, . . . , N , and separated from the remainder of the spectrum 
(such a contour exists in view of Assumption [2]) , and define 

FN+i{t) = -^l {y,Pk{t)\ R{z,t)R{z,t)dz, (3.23) 
GN+i{t) = -^i i?(z,t)i2(z,t) ( VPfc(t) ) dz. (3.24) 



By definition of the contours. 



--^ / R{z,t)dz = Pj{t), l<j<N, 



and 

N 



~ I Riz,t)dz = Y,Pk{t) = P^+i{t). 



Besides, in view of the continuity of 1 1-^ Ej{t) for all 1 < j < A^, it is possible to use contours 
which are locally constant in time, i. e. for a given t > — oo such that f{t) > 0, there exists a 
(small) time interval {t — T,t + t) and a contour F*- such that 

Vs G (t - T,t + r), — (h R{z,s)dz = Pj{s) 

2ivr Jr' 

for 1 < i < A^, a similar result holding for j = A^ + 1. Using such locally constant contours, 
the time derivative of the contour integral defining the projector can be restated as a contour 
integral of the time derivative of the resolvent: 

_1 / R^,^t)dz = ^, 1<J<N, 
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and 



— <f> R(z,t)dz = V , , . 

Boundedness of F. The operator F{t) is bounded. To see this, we first rewrite Fj (1 < 

j < N) as 



Fj{t) = P^{t)R{E,{t),t)Pf{t)^^. (3.25) 
Indeed, using the expression (|3.2ip of Fj, 

F,{t)-PHt)R{E,{t),t)PHt)'^^'^^^ 



^ I P^^{t){R{z,t) - R{Ej{t),t))P^^{t)R{z,t)dz 

— / PHt){R{z,t) - R{EAt),t))PHt)R{z,t)H{t)R{z,t)dz. 
2ivr JvAt') 



lT,{t) 

i 

When the contour encircles closely enough Ej{t) 



\\Riz,t)\\< ^ 



6, it) 



Using the resolvent identity, it follows 



\P^^it)iR{z,t) - RiE,it),t))Pf{t)R{z,t)\\ = \z-Ej{t)\ ■ \\R{z,t)Pf{t)R{Ej{t),t)P^{t)R{z, 

\z-Ej{t)\ 



< 



Then, the difference 

/ PHt){R{z,t) - R{EAt),t))R{z,t)PHt)H{t)R{z,t)dz 



fit) 



can be made arbitrarily small by decreasing the radius of the contour Tj{t), with a constant 
C depending on the relative i?o-bound of V. 

From the expression ([S^SI), and the bound \\Pj^{t)R{Ej{t),t)P^^{t)\\ < 5j{t)-^, it holds 
finally 

iif,mii<iaM<c^'« 



where we recall that both / and /' are non-negative. This shows that Fj{t) is a bounded 
operator when f{t) > 0. A similar bound holds for Gj. 

The terms F/v_|_i(t), require a different treatment. In this case, the uniformity of 

the gap between the N eigenvalues encircled by Tj\f^i{t), and the remainder of the spectrum 
may be used to construct a contour T]\f^i{t) such that 

yz€TN+iit), \\R{z,t)\\< ^ 



A{ty 
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This can be done by ensuring that the contour remains far away enough from the remainder of 
the spectrum, while still being at a finite distance of the first N eigenvalues. In particular, it is 
possible to construct a contour intersecting the real axis at a point 7 such that I7 — Eiy{t)\ > 
A(t)/4 and 

inf{|7-E|, E£a{H{t))\{Ei{t),...,EN{t)}}>A{t)/A. 

Then, 



\\Fn+i 



fit) 



2i7r 



N 



VPfc(t) R{z,tfVR{z,t)dz 



< c 



A(t)3' 



(3.26) 



and so F/v+i is bounded since A(t) > A* > and /' is bounded. A similar bound holds for 

Gn+i- 

In conclusion, 

<Cf^w, (3.27) 



fit)' 



for some constant Cp independent of t. 



Computation of the commutator. It is easily seen that F{t) maps the Hilbert space TC 
to D{Hq). The commutator [H{t), F(t)] can then be defined as an unbounded operator with 
domain D{H{t)) = D{Hq). For a given 1 < j < A^, it holds, using the commutation property 



mt),F,{t)] 







2i7r 


/ 

Jrjit) 






2i7r 








2i7r 


/ 



[H{t),PHt)R{z,t)R{z,t)]dz 



[H{t) - z, P^{t)R{z, t)R{z, t)] dz 



-PHt) 



dPj{t) 



dt 



2m 



+ ^^WU-9 R{z,tfdz\H{t) 



v,{t) 



-{l-P,{t)) 



dt 



following the proof of Theorem 2.2 in [22]. Similar computations show 



mt),Gj{t)] 



dPj{t) 
dt 



(I-P,(t)). 



Finally, for 1 < j < iV, 

[H{t),F,{t) + G,{t)] = 

In the same way, 

[H{t),FN+l{t) + GN+l{t)] = - PN+l{t), 



Pjit), 



dPj{t) 



dt 



dt 



PN+i{t), 



dPN+ijt) 
dt 
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Since 



Kit) 



N+l 



dPjit) 
dt 



-I N+l 

-Y 



dt 



it holds 



[if(t),F(t)]=i^(t). 



(3.28) 



Integration by parts. Consider now — oo < to < ^ ^ such that /(to) > (hence f{t) > 
since / is non-decreasing). Define 

Kit) = -ieU%to,t)F{t)U'{t,to). 



Then 



JC'it) = U%to,t)[Hit),Fit)]W{t,to)-ieW{to,t)F'it)U%t,to). 



In view of (I3.28p . the difference between the evolution operators is rewritten as 



U'{t,to)-Ui{t,to) = -U'{t,to) U'{to,t')Kit')Uiit',to)dt' 

J to 



-U'HM (^^^ + ieU^to,t')^^U%t',to)^ U%to,t')Ui{t',to)dt', 

(3.29) 



so that, after an integration by parts for the term associated with /C', 



\U'{t,to)-Ulit,to)\\ 



U'{to,t')Kit')Ulit',to)dt' 



to 

< ||/C(t)|| + ||/C(to)||+e /V' 

J to 

rt 



(3.30) 



' nt')^{u'M)mt'M) dt' 



to 



dt' 

t 



< e(\\F{t)\\ + \\F{to)\\+ [ \\F'it')\\dt'+ [ \\F{t')\\\\Kit')\\dt').i3.3l) 

\ Jto Jto / 



The first two terms in the above equality are bounded with the bound (|3.27|) on F. For the 
last one, we use again the bound (|3.27|) on F, and the fact that K is uniformly bounded 
(see dSJlD), so that 



\\F{t')\\\\K{t')\\dt' <C f < ^ [\f'f. 
to Jto J JK'^O) Jto 



(3.32) 
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We now turn to the central term. For 1 < J < A^, and using (|3.25|) . 



\F'{t')\\dt' < 



to 



< 



< 



< 



*^<^<'./'ll^(.')ll ' 

to "j \^ ) J to 



to 



to 



dt 
dt' 



-[P^{t')R{E,{t'),t')P^{t') 



dt' 



+ 



to 



to 



< c 



to 



fit') 



fit') 



W) ' \f{t')) 



for some constants C, C" > (related to the relative ffo-bound of V). Similar expressions 
can be obtained for Gj (1 < j < N). Straightforward estimates can be used for F/v+i, Gn+i-, 
following a treatment similar to what was done to obtain (|3.26|) . upon deriving the terms 
appearing in the contour integral: 



with 

In conclusion, 



f'{t) dxP{f{t)) <Gf'{t) 



for some constant C > 0. 



(3.33) 



Decomposition of the integral close to the degeneracy. In order to avoid the singulari- 
ties when /(to) 0, the difference of the unitary operators is separated into two contributions 
as 



[/^(0,to)-C/i(0,to) 



-C/^(0,to) / U'ito,t)K{t)UUt,to)dt 

J to 

-U%0, to) U'{to,t)K{t)Ul{t, to) dt, 

where T is chosen such that f{T) > 0. The first term is bounded using the straightforward 
estimate 



U'iO,to) I U%to,t)K{t)Ulit,to)dt 

to 



<C [ Y.mmdt<C' [ f'{t)dt<G'f{T). 

Jto Jto 

(3.34) 
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For t E [T, 0], f{t) > f{T) > and there is a gap proportional to f{T) between the eigenvalues: 

Vl<i<iV, VtG[0,T], 6,{t)>af{T), 



for some a > 0. The inequality (|3.3Q|) . combined with (|3.27p . (|3.32|) and (|3.33p . allows to 
bound the second term as 



U'{0,to) / U%to,t)Kit)Ulit,to)dt 

JT 

- V/(o) f{T) 



T 



U'{T, t)K{t)Ul{t,T)dt 



(3.35) 



The limit to ^ — oo can then be taken in the above expressions. Moreover, upon choosing T 
small enough so that f{T) = 6^/3 < 1, it follows, adding (pOll and jS^Sl), and using the fact 
that /' G Li((-oo,0]) nL°°((-oo,0]) and /" G Li((-oo,0]), 



|[/-(0,-oo)-C/i(0,-oo)|| <C /(r) + e 1 + 



(3.36) 



This concludes the proof. 



3.4 Extensions 



The above proofs can be straightforwardly extended to the following cases (see Section [2] for 
the notation). 



Definition of the initial states when VqVVq has degenerate eigenvalues. Two changes 
should be made in the proofs presented in this paper: (i) the estimate obtained in the adiabatic 
limit degrades; (ii) more conditions are required to define the initial states. 

Denote by £"0,1 the M < N eigenspaces associated with the eigenvalues of VqWo, set 
rii = dim(£'o,i), and define 



A/; = j/c G {!,..., iV} (/?fc,o e fo,i} 



the set of indices corresponding to the i-th eigenspace of VoVVq. Of course, 

M 



^ni = iV, Card(7Vi) = nj. 



In view of Assumption [3l for any {k,l) G J\f^, k ^ I, there exists an integer pk^i > 2 and an 
analytic function eki{\) such that 

Ek{X) - Ei{X) = XP^'^ CkiiX), efc,z(0) / 0. 



Denote by p^, the maximal integer for all couples 1 < k,l < N . Then, the final estimate (|3.36p 
in the proof of the adiabatic limit reads 



\U^0,-oo) - C/i(0,-oo)|| < 



c(/(r) + .(i + -^))<3CeV(2..+i), 
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which is indeed larger than the e^/^ bound found in the case p = 1 (no degeneracy of the 
perturbation restricted to Sq). 

We now describe an iterative procedure which determines the initial states in a unique 
manner, using the higher order equations in the hierarchy (|3.10p . We start with the conditions 
of order 2. A necessary condition for (|3.1Q|) to have a solution is that its right-hand side belongs 
to £"0^. With (Cnal) . this requires, for all l<j,k<N, 

{ ipk,o, VRoVipjfl ) + Ej^26j,k + {Ej,i - Ek,i)c]^f, = 0, (3.37) 

where J^^b is the Kronecker symbol. In particular, 

Vi €{!,..., M}, y{j,k)eAf,^, {ipk,o,VRoVipj,o) + Ej,2Sj,k = 0- 

Therefore, {(/^j^oIjgM has to be an eigenbasis of Vo^iVRoWo^i where Vo^i denotes the pro- 
jector onto (?o,i- If 'Po,iy RoVPo^i has non-degenerate eigenvalues, the initial eigenfunctions 
{fk,o}k£j\fi are uniquely defined. 

Otherwise, the procedure must be repeated. Recall that there exists an integer such 
that after steps the degeneracy has no further split (see the discussion at the beginning of 
this paragraph). When the degeneracy is not permanent (see below for this case), this allows 
to determine the initial states in a unique manner. See for instance In many practical 
cases however, degeneracy is never totally split because V shares some symmetries with Hq. 
In this case, permanent degeneracy has to be taken into account (see below). 

Decomposition of the switching. In the case when (|2.2|) is not satisfied, i.e. ||-P;(0) — 
-P(— oo)|| = 1 or equivalently ||Pj(0)^/^j || = (since the eigenspaces are assumed to be one- 
dimensional), the switching should be done in several steps. The intermediate steps can be 
chosen by finding a finite number of values G [0, 1] (A; = 1, . . . , — 1), with Aq = and 
Aat = 1, such that ||Pj(Afc+i) — Pj{Xk)\\ < 1. This is possible since Pj is continuous on the 
compact interval [0, 1] . 

The initial state ipo is evolved into a state ^pl by switching from Hq to Hq + XiV as 

ipi = hm —. — 
where the evolution operator 



is the following operator in the interaction picture: 

dm (t,to) / \ 

ie^^^ = [Ho + Xif{t)VjUl^{t,to), C/|^(to,to) = I- 

The state ^pl is then evolved into a state tp2 by switching Hq + XiV to Hq + X2V as 



ip2 = lim ■ 
where the evolution operator 
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is defined as the following operator in the interaction picture: 

ie '^^'^'^f = {Ho + X,V + {X2-Xi)fit)v)ui^^^^{t,to), ^A„A,(io, to) = L 

This construction is repeated until an eigenstate ipN of i^o + ^ = Hq + Xj\fV is obtained. 
Notice that it is important to do the procedure sequentially. 

Permanently degenerate eigenspaces. When there are permanently degenerate eigenspaces 
associated with one of the eigenvalues Ej{X) or Ej{t), the determination of the initial basis can 
still be performed as it is presented in Section 13.21 However, the argument leading to (I3.19P 
in Section [3.3. II cannot be extended as such to the case when RanPj (O) is of dimension larger 
or equal to 2. This is not a problem since ^4(0, —oo)ipj is still an eigenvector of Pj{0), and its 
phase can be removed upon considering 

f^l,int(0'-H^i _ ^(0,-00)^J 



i<t>\ U%,,JO,-oo)i;j) {^\ A{0,-oo)i;j) 

for some fixed state (/>, provided the denominator is non zero. In Theorem [TJ the choice (j) = V'j 
is done, together with the assumption ( cj) \ ^(0, — oo)^j ) / 0. This assumption could in this 
specific case be translated into an assumption on ||-Pj(0) — Pj(— oo)||, but in general it should 
then be assumed that there exists cp £ 7i such that ( (p \ A{0, —oo)tpj ) 7^ 0. 



Existence of finitely many crossings. The projectors being analytic, the Kato operator 
can still be defined when there are eigenvalue crossings. The main issue in extending the 
Gell-Mann and Low formula to this case is therefore the proof of the adiabatic limit, which 
can however still be handled with |22l Corollary 2.5] since the crossings are regular (again, 
because the eigenvalues are analytic). 



Initial subspace composed of several degenerate spaces Sq, £1,... In this case, the 
operator V should be diagonalized in each subspaces, i. e. the self-adjoint finite-rank operators 
VjVVj\^ are diagonalized in order to construct a basis of £j. A global basis is then obtained 
by concatenation (direct sum). 
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